Intrinsic and extrinsic mirror symmetry breaking in anti-dot spin-wave waveguides 
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We theoretically study the spin-wave spectra in magnonic waveguides periodically patterned with 
square anti-dots in nanoscale with pinned magnetization at the edges. We show that the breaking 
of the mirror symmetry of the waveguide by the structural changes can result in a magnonic band 
gap closing. These intrinsic symmetry breaking can be compensated by properly chosen asymmetric 
external bias magnetic field, i.e. in an extrinsic way. As a result the magnonic gaps existing in the 
ideal symmetric structure can be recovered. The model used for the explanation also suggests that 
this idea could be generic both for exchange and dipolar interaction regimes of spin-waves and also 
for other types of waves, e.g., electrons in the graphene ribbons. 

PACS numbers: 75.30.Ds,75.75.-c,75.78.Cd,75.70.Ak 



The possibilities for fabrication metallic magnetic ma- 
terials with nanoscale precision has opened the way for 
tailoring the dispersion of high-frequency spin waves 
(SWs) . This can be done also with the use of magnonic 
crystals (MCs), which are magnetic counterparts of pho- 
tonic crystals. [TJ [2] Recently the first bi-component MCs 
have been realized at the nanoscale and the opening of 
magnonic gaps were proved experimentally. [51 0] Two- 
dimensional (2D) anti-dots lattices (ADL) formed by pe- 
riodic array of holes in ferromagnetic film can be much 
easily fabricated. These systems have been intensively 
studied in recent years on length scales from microm- 
eters down to few tens of nanometers. [5HT0] For ADLs 
with big period the inhomogeneity of the internal mag- 
netic field is decisive for formation of the magnonic band 
structure. [TT] With the decreasing period of ADL the 
Brillouin zone (BZ) border will move to larger wave- 
vectors and the exchange interactions at some point start 
to play a primary role in the formation of a magnonic 
band structure. 

Prototypes of basic magnonic devices have already 
been demonstrated to be promising for technological 
applications [T2TH3] but the scaling down of magnonic el- 
ements to tens of nanometers in size and tens of GHz op- 
erating frequencies [T5T - H7] are still a challenge. Waveg- 
uides for SWs, which are an important component in 
most magnonic devices. |12j have been realized experi- 
mentally so far only for SWs in the frequency range up 
to few GHz.[7l lT8H2"2"] To predict properties of magnonic 
devices at scale of tens of nm, more basic research needs 
to be conducted. Therefore, theoretical investigation of 
the SW waveguides operating in the range of tens of GHz 
is a frontier field of research. 

Recently the micromagnetic simulations (MSs) were 
used to show that periodic waveguides have filter prop- 
erties due to the opening of magnonic gaps in the SW 
spectrum. 23 26] In this Letter we investigate how loss 
of a mirror symmetry within a one-dimensional (ID) 



nanoscale anti-dot lattice waveguide (ADLW) may affect 
the SW dispersion. When this symmetry exists, based 
on their profile with respect to the central longitudinal 
axis, the SW spectra can be separated into two groups: 
symmetric modes and anti-symmetric modes. The break- 
ing of the mirror symmetry will automatically make the 
classification impossible. We will study two types of the 
symmetry breaking mechanisms namely, intrinsic and ex- 
trinsic, which are realized by changes of the position of 
the anti-dots row and by application of asymmetric bias 
magnetic field, respectively. The question is how do these 
changes influence the magnonic spectra and the existing 
magnonic gaps? How big the symmetry breaking needs 
to be in order to close the gap? We perform theoret- 
ical calculations with the plane wave method (PWM) 
and confirm the results by using MSs, to address these 
questions. The discussion on the possible extensions of 
our findings to large scale and other systems will be con- 
ducted on basis of an analytical model. 

We study the symmetric and asymmetric magnonic 
waveguides based on the anti-dots lattice structure shown 
in Fig.nl It has a form of a thin (thickness u — 3 nm) and 
infinitely long permalloy (Nig Fe2o) stripe with a single 
row of square holes (s = 6 nm anti-dots) disposed peri- 
odically along the central line. The stripe width and the 
lattice constant are fixed at 2 x w + s = 45 nm and a = 
15 nm, respectively. The row of holes is placed in the 
distance w — 19.5 nm from the top (and bottom) edge 
of the stripe in the case of the symmetric ADLW. A bias 
magnetic field is applied along the stripe and it is strong 
enough to saturate the sample (/xo-ffo = 1 T). The satu- 
ration magnetization M s = 0.8 x 10 6 A/m, the exchange 
length A cx = 5.69 nm, and the gyromagnetic ratio 7 = 
175.9 GHz/T were assumed in calculations. 

The calculations of the magnonic band structure are 
performed with the finite difference method MS and 
the PWM, with OOMMF[27] and a Fortran code devel- 
oped by authors, respectively. Both methods solve the 




FIG. 1: Anti-dots lattice waveguide under investigation: u = 
3 nm thick and 2w + s = 45 nm wide infinite Py stripe with 
a periodic series of square anti-dots (of size s x s, where s = 
6 nm) disposed along the waveguide with a period of a — 15 
nm. The row of anti-dots divides the waveguide into two sub- 
waveguides of the width w = 19.5 nm each. Bias magnetic 
field fioHo = 1 T is oriented along the waveguide. 
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Landau-Lifshitz-Gilbcrt (LLG) equation: 
dM(r,t) 



m 



7jU M(r,t) xH eff (r,t) 
<*{■*., ■< dM(r,t) 



(1) 



where r and t are position vector and time, respectively. 
Ho is the permeability of vacuum. The first term on the 
right hand side is related to the torque inducing preces- 
sion of the magnetization M and the second one describe 
the damping process (a denotes the damping constant). 
The damping is neglected in PWM calculations and in- 
cluded in MS (a = 0.0001). The effective magnetic field 
H e fj here consists of the bias magnetic field Ho, demag- 
netizing field and exchange field. The pinned dynamical 
components of the magnetization vector were assumed at 
Py/air interfaces in calculations with both methods. The 
pinning in OOMMF was introduced by fixing magneti- 
zation vector in all cells of the discretization mesh which 
border the anti-dots, i.e., for the width 0.5 nm along y 
axis. (In MS the discrete mesh size of 1.5 x 0.5 x 3 nm 
along x, y and z axis, respectively, were used. The MS 
were performed for 4 ns. In the PWM we use 961 plane 
waves.) In the PWM the pinning is applied exactly at 
the edges of Py. Due to small thickness of the ADLW, 
uniform SW profile across the thickness is assumed. Both 
methods were already used in the calculations of the SW 
dynamics and proved to give correct results. [51 [HI [2S] 

We start our investigation for the symmetric ADLW 
(Fig. |T|).[28 The dispersion relation of SWs in the sym- 
metric ADLW is shown in Fig. [2] the results of the 
OOMMF simulations are shown in Fig. W[&) and of the 
PWM in Fig.[2Fc) (black solid lines). The agreement be- 
tween results of these two methods is satisfactory. The 
presence of two magnonic band gaps is evident and they 
are marked by the yellow color. The origin of these two 
band gaps was found to be different. The first one opens 
at the BZ boundary due to a Bragg reflection of SWs, 
while the second gap is opened in the middle of the 



FIG. 2: (Color online) Magnonic band structure of ADLWs 
(shown in insets above the main figures where the thin dashed 
line marks the middle of the ADLW) calculated with OOMMF 
in (a) and (b), and with PWM in (c). The band structure for 
the symmetric ADLW are shown in (a) and in (c) with black 
solid lines. The results for ADLW with shifted upward anti- 
dots row are shown in (b) and in (c) with red lines calculated 
with OOMMF (Aw = 1 nm) and PWM (Aw = 0.9 nm), 
respectively. At the bottom the squared amplitudes of the 
dynamical magnetization |m z | 2 for first four modes in the 
center of the BZ (calculated with PWM-cf. (c)) are shown 
for symmetric and asymmetric ADLW, first and the second 
row, respectively. 



BZ.J2SJ It was shown that this splitting of the bands is 
due to the anti-crossing between two family of modes, |23j 
those with and without a nodal line in the upper and 
lower part of the ADLW (see the first row of profiles in 
Fig.§. 

The structure investigated above has a mirror symme- 
try with respect to the central axis of the ADLW. The 
1 st and the 2 nd mode (and also 3 rd and 4 th ) is symmet- 
ric and antisymmetric with respect to the central axis of 
the ADLW, respectively (see the first row of profiles at 
the bottom in Fig. [2]) and their frequencies are degener- 
ate. The shift of the row of anti-dots from the central line 
will break this symmetry. Frequencies of modes will split, 
with one mode shifted up and the second mode shifted 
down in the frequency scale. The dispersion relations of 
SWs in asymmetric ADLWs, obtained by shifting the row 
of anti-dots by Aw = 1 nm and 0.9 nm upward and cal- 
culated using OOMMF and PWM, are respectively pre- 
sented in Fig. [2jb) and in (c) (red lines). (In OOMMF 
slightly larger value of Aw were used because of the lim- 
itations of the discretization mesh and time needed for 
simulations.) We see that the shift of the anti-dots row 
by only 2% of the total width of the ADLW is enough 
to close both magnonic gaps. The first mode, with lower 
frequency, is now concentrated in the wider part of the 
ADLW (width w + Aw) (see, Fig. [2] in the bottom). The 



second mode, with higher frequency, has an amplitude 
concentreted in the upper part of the ADLW (of width 
w — Aw). The same mechanism is responsible also for 
closing the second gap, even though the origin of this gap 
is different and the respective changes of the frequencies 
of the third and the fourth bands are larger. 
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FIG. 3: (Color online) Width of magnonic gaps in ADLW 
as a function of (a) the shift of the anti-dots row for Aid 
and (b) an additional asymmetric bias magnetic field A7?o 
in the symmetric ADLW. The AHq increases bias magnetic 
field in the upper half of ADLW (to the value Ho + AHq) and 
decreases bias magnetic field below the row of anti-dots (to 
the value Ho — A Ho). 

The gap widths as a function of the value of Aw are 
shown in Fig. f3^a). The widths decrease monotonously 
with increasing Aw, the slope for the second gap is bigger 
and this gap is closed at 0.45 nm, while the first gap exists 
up to 0.8 nm. 

We have just shown that the magnonic spectra and 
especially the magnonic gaps in the considered ADLW 
are very sensitive to a mirror symmetry breaking. The 
question comes: is it possible to compensate the effect of 
an intrinsic symmetry breaking in the ADLW by the ex- 
trinsic one? In our case it will be a compensation of the 
effect of the structural asymmetry on the magnonic band 
structure (and magnonic gaps) by asymmetric bias mag- 
netic field. The answer will begin from the consideration 
of the analytical model. 

In the homogeneous waveguide (WG) the solutions of 
the linearized LLG equation (with damping neglected) 
can be written in the following form: m(x,y) — 
m(y)e xX , where k x is the wave- vector of the SW along 
the WG and m(y) describe the dependence of the am- 
plitude of dynamical components of the magnetization 
m across the WG width (we assume the uniform mag- 
netization across the WG thickness which is much less 
than the width). The solutions can be estimated as: 
m(y) f=a sin(«:j/), cos(Ky) where the transversal compo- 
nent of the wave-vector k, — (n + l)iv/w e g is quantized 
{n = 0, 1, 2, ... counts the number of nodal lines across 
the WG width). For strong but not ideal pinning the 



effective width w e g — wd/(d — 2) depends on the pin- 
ning parameter d, which determines the boundary con- 
ditions (BC) for magnetization and gives also a possi- 
bility to include the dipolar effects into the model. [3TJ] 
It varies in general from to oo for the transition 
from unpinned to pinned BC. The pinning parameter 
d = ^(l-Jj^)/m~2M%) + (^) 2 )} depends both 
on the material and structural parameters (K s denotes 
the surface anisotropy). For d w oo, as in our numerical 
calculations, n = means no nodal line in the upper or 
lower part of ADLW (see, Fig. pi modes 1 and 2), n = I 
a single nodal line (see, Figpj modes 3 and 4), etc. The 
dispersion relation of SWs in the WG can be written in 
the form:[ 



LU = \J (Wo + W cx ) (Wo + "ex + Wdip), 

where lu is an angular frequency of SWs. u>q = ^jjiqHq, 

u cx = 7MoM s |f (nV + k 2 x w 2 ), w dip = 7 /; Af s 1 ~ c % ( u fc '" ) 
denote the contributions from external, exchange and 
dipolar fields, respectively. 

The estimations of changes in SW dispersion relation 
resulting from the changes of w or H$ can be done by 
calculation of the full differential of the function oj = 
u>(w,Hq). It allows to derive the relation between small 
changes of AHq and Aw, for which the compensation of 
intrinsic and extrinsic symmetry breaking can occur, i.e., 
when duj(w, Hq) = 0: 



HoAHo 2n 2 fj, M s X 2 x {n + l) 



Aw 



/' 



K a 



irufioM 2 



(2) 



This ratio has units of T/m and shows how big asym- 
metric magnetic field needs to be added to compensate 
the shift of the anti-dots. The function: f(s,l,r) — 
[s + I(2£ - 4rln(r))][s + I(r + £ - 2rln(r)) - S )]s~ 2 
depends on three dimensionless parameters: s-the rela- 
tive strength of the surface anisotropy, ^-ratio between 
exchange length and the width of WG and r-the aspect 
ratio of the WG. The values of f(s, I, r) with big absolute 
value of K s , refer to the regime of strong pinning. Note 
that Eq. (pi) does not depend on k x , which means that it 
should be fulfilled for any wave-vector. 

In our ADLW, we have two WGs separated by the anti- 
dots row. When we shift the row of the anti-dots Aw 
along positive y direction, the width of the upper WG is 
decreasing by Aw and the width of the lower WG is in- 
creasing by the same amount. It will result in an increase 
of the frequency of the modes concentrated in the upper 
WG and decrease of the frequency of the mode concen- 
trated in the lower WG. To compensate for these changes 
in the dispersion relation by a bias magnetic field we need 
to do the opposite. According to Eq. (pi) we need to ap- 
ply different bias magnetic fields to upper and lower WG. 



The dependence of the magnonic gap width under asym- 
metric bias magnetic field applied to the ADLW, AHq 
(i.e., in the upper part of the ADLW bias magnetic field is 
H + AH , while in the lower part of ADLW is H -AH ) 
calculated with PWM is shown in the Fig. [3Tb). 

After these estimations we perform the PWM calcula- 
tions. The results are presented in the Fig. Bb) and (c) 
for AH to recover the first and the second magnonic gap 
in the asymmetric ADLW (i.e., when Aw — 0.9 nm and 
0.5 nm), respectively (see FigJ3^i). It is interesting that it 
is possible to recover the first and the second magnonic 
gaps but with different values of the ratio Mo Att) " . The 
analytical values of this ratio (calculated from the Eq. (pi) 
with w = 18.5 nm, i.e., the distance between pinned lay- 
ers used in MS) for the ideal pinning (f(s,l,r) = 1) for 
the first gap (when n = 0) and the second one (n = 1) are 
101 mT/nm and 406 mT/nm, respectively. To validate 
our predictions we performed MSs for /iqAHq = 105 mT 
and Aw = 1 nm. The MS results are shown in Fig. (4th) 
with the first frequency gap opened and with good agree- 
ment with PWM calculations shown in Fig. Efb). In or- 
der to open the second gap we have to cancel the shifts 
for the bands with a single nodal line (n = 1) (this gap 
is formed due to the anti-crossing of the n = and n = 1 
modes but at the Aw = 0.5 nm the splitting of the n = 1 
dominate) by applying the field for which ^ Aw ° is about 
4 times bigger than for the first gap (410 mT/nm). This 
confirms the applicability of Eq. (pi) with a square de- 
pendence on n + 1. The profiles of SWs (compare bot- 
tom panels of Fig. p] and Fig. |4| further establish the 
restoration of amplitude distribution by extrinsic com- 
pensation. Small differences in the extent of band gap 
recovery obtained from numerical calculations and the 
analytical model, show that the pinning in the middle of 
ADLW is not perfect. 

We have checked also that the compensation effect 
will still appear for continuous change of the bias mag- 
netic field, which is also more justified in an experimen- 
tal research. We used ramp-like profile of the magnetic 
field: [32] H = H Q + 2AH Q {2y + Aw) where y = corre- 
sponds to the ADLW center. The values of Aw and AHq 
are the same as for step-like profile of magnetic field. We 
have found, that this kind of the field acts similarly to the 
step-like profile of the field, when its value is normalize 
to have the same average value for corresponding sub- 
waveguides (the formula for ramp-like field fulfills this 
condition) . The results of PWM calculations are shown 
in Fig. [4] in right insets, this was also confirmed in MSs. 

We have shown that a small symmetry braking in anti- 
dots waveguide can result in closing of magnonic gaps in 
the range of the spectra determined mainly by exchange 
interactions. It can be crucial for practical realization 
of the SW waveguides because precise mirror symmetry 
will be difficult to achieve on such small scales. We have 
shown that the effect of a deviation from the ideal struc- 
ture in magnonic band structure can be compensated by 
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FIG. 4: (Color online) Magnonic band structure of the ADLW 
showing the compensation effect of the intrinsic asymmetry by 
the extrinsic field. In (a) and (b) (green dashed lines) the first 
gap (for the modes without a nodal line — n — 0) is completely 
reopened. The calculations with OOMMF (a) and PWM (b) 
were performed for Aw = 1 nm and 0.9 nm, respectively, 
with fj, AH = 105 mT (fi H = 1 T). The reopening of the 
second gap (opened in the anti-crossing of the mode n = 
with n = 1) is presented in (c). Calculations in (c) were done 
with PWM for fi AH = 205 mT and Aw = 0.5 nm. The 
left insets in (b) and (c) show enlarged results for the step- 
like field profile of the bias magnetic field, the right ones - the 
outcomes for linear change of the magnetic field profile (ramp- 
like profile) across the ADWL. At the bottom row, profiles of 
SW calculated with PWM are shown. Profiles for modes 1 
and 2 are calculated for the structure from (b) and modes 3 
and 4 for the structure from (c) at the BZ center. 



the bias magnetic field, i.e., in an extrinsic way. Using the 
predictions of an analytical model, we showed in numeri- 
cal calculations that such a compensation is possible but 
the field should be asymmetric as well. It was presented 
that two magnonic band gaps of different origins can be 
reopened in asymmetric waveguides by this way. Our 
predictions are based on the symmetry properties and a 
general analytical model, so it should be possible to ex- 
tend this idea to larger structures, i.e., closer to actual 
experimental samples where the dipolar interactions are 
important. The intrinsic and extrinsic symmetry break- 
ing or its compensation can be also exploited to manip- 
ulate active and inactive modes, which couple to the ex- 
ternal fields, [33, 34 in a similar way as was predicted for 
plasmonic metamaterials. [35] Further study is necessary 
to investigate the compensation effects proposed here in 
other systems, e.g., like in electrons propagating in an 
periodically patterned graphene ribbons by the external 
electric field. [31 [37] 
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